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Random banded matrices with linearly increasing diago- 
nal elements are recently considered as an attractive model 
for complex nuclei and atoms. Apart from early papers by 
Wigner there were no analytical studies on the subject. 
In this letter we present analytical and numerical results for 
local spectral density of states (LDOS) for more general case 
of matrices with a sparsity inside the band. The crossover 
from the semicircle form of LDOS to that given by the Breit- 
Wigner formula is studied in detail. 

PACS numbers: 05.45. +b 



Recently there was a growing interest in the statisti- 
] cal properties of Random Banded Matrices (RBM) (see, 
e.g. and references therein). These N x N matrices 
Hnm can be characterized as those containing nonzero el- 
ements only within a wide band of the size b ^ N around 
the main diagonal. It was realized that matrices with 
such a structure are appreciably different in their statis- 
tical properties from those forming the classical Gaussian 
' Ensembles studied by many authors j^. More precisely, 
the RBM of infinite size do not show the famous effect 
of level repulsion that is the direct consequence of the 
localized nature of their eigenvectors . 

The most studied type of RBM is that with the zero 
mean value of all matrix elements and with the variance 
' given by i/,^„„ — V^ad n — m \ /b), the function a(r) 
decaying exponentially (or faster) at r 3> 1. Such ma- 
trices were shown to be relevant for understanding some 
properties of the Quantum Kicked Rotator, one of the 
aradigmatic models in the domain of Quantum Chaos 
On the other hand, the matrix Hnm can be treated 
both as a tight-binding Hamiltonian of a quantum parti- 
cle in Id system with long range random hoppings and as 
an adequate model for quasi Id disordered wires 0|. In a 
number of publications properties of such matrices were 
studied both numerically and analytically, see references 
in 1^. In particular, it was found that for matrices of in- 
finite size any eigenvector has typically finite number (of 
the order of /oo oc &^) of essentially nonzero components. 
When the matrix is of a large finite size iV ^ 6 ^ 1, its 
statistical properties were shown to be determined by the 



scaling parameter b'^/N As to the density of states 

(DOS) it has generally a form of the semicircular law [pj: 
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with i?o oc 6^/^ being the semicircle radius. 

Another type of the RBM - that with the mean value 
of diagonal elements increasing linearly along the main 
diagonal: 77„„ = an - has attracted recently a lot of 
research activity . In a tight-binding analogy men- 

tioned above these matrices describe a quantum particle 
in Id disordered system subject to a constant electric 
field. Another interpretation was developed in [pO 12 
where matrices of such kind are considered as Hamilto- 
nians of generic conservative systems with complex be- 
haviour, like heavy atoms and nuclei. Indeed, a very de- 
tailed study of compound states in the chaotic Ce atom 
Jl^ , p^ has revealed that the latter class of the RBM can 
be used rather efficiently to describe such a physical sys- 
tem. This particular kind of RBM was introduced for 
the first time by Wigner about forty years ago ||l|. For 
this reason we refer to this ensemble as that formed by 
Wigner Random Banded Matrices (WRBM). Quite close 
ideas were developed recently ||T^ where the model for 
a general integrable Hamiltonian Hq perturbed by some 
generic perturbation Vo was investigated. It was argued 
that the matrix {Hq + VQ)nm in the basis of semiclassical 
eigenstates of the unperturbed part has the form of 
sparse WRBM. 

As is known, much more informative characteristics 
compared to the global DOS is the so-called local spec- 
tral density oj states (in the literature, LDOS). The latter 
function, also known in nuclear physics as strength func- 
tion, is defined as 



p(£;,n)=^ |7^,(n) ^'SiE-E, 
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where ^pv{n) is n-th component of the eigenvector ipi, 
corresponding to the energy level E^, n = 1,2, ...N. Af- 
ter averaging over the disorder the LDOS depends (in 
the limit N — > oo) on the parameter z = E — an only: 
p{E,n) = pl{z). As a result, the global DOS is energy- 
independent: p{E) = provided 
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the parameter a is small enough for the sum being re- 
placed by the integral |p^ . 

The quantity of the most physical interest is the width 
r of the LDOS: it generates a new length scale Imax = 
T/a, see details in 10 1^. If one treats linearly increas- 
ing diagonal elements as eigenvalues of an unperturbed 
Hamiltonian the off-diagonal elements play a role of per- 
turbation and the length Imax can be interpreted as an 
effective number of unperturbed states coupled by the 
perturbation |p^ . In the limit Imax 1 standard per- 
turbation theory can be applied since the perturbation is 
weak. We will refer to this regime as the perturbative one. 
In the opposite limit Imax S> 1 the problem is essentially 
non-perturbative. The latter case is our main concern in 
the present paper. 

In the case of very large a well inside the perturbative 
regime the form of LDOS is determined mostly by the 
distribution of diagonal elements (e.g. the Gaussian). 
On the other hand, for a in the nonperturbative regime 
the region of parameters was indicated in |17| where the 
Lorentzian shape 



PBwiE,n) = 
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is expected with F « 2TrpV^ jl^. Actually, this form 
of LDOS was obtained by Wigner in his early studies of 
WRBM as a result of some limiting procedure. 

The equation Eq.(||) is commonly accepted approxi- 
mation for the LDOS in nuclear physics known as the 
Breit- Wigner (BW) formula. It is considered to be well 
in agreement with experimental data for nuclei and 
complex atoms [Q. This fact makes the WRBM to be 
much more attractive as a model for physical systems 
as compared with the full random matrices. However, a 
broader application of this ensemble is restricted by the 
absence of a general theoretical understanding of their 
properties. Indeed, the method used in Wigner pioneer- 
ing papers |^ is rather involved, partially heuristic and 
seems to be based upon a specific form of the parameters 
chosen. In particular, it is quite unclear how universal 
are Wigner results (e.g., the Breit- Wigner form of the 
LDOS, Eq.(^) against variations of the parameters a, b, 
form of the distribution of the off-diagonal elements, and, 
ultimately, the sparsity po| . 

Motivated by all this, we have performed an analytical 
consideration of the LDOS for general WRBM with ar- 
bitrary degree of sparsity. For this purpose we consider 
random symmetric matrices Hnm = OLnSnm+inm , where 
the probability distribution of , n < m is given by: 



nm) H~ J^nm^l (^nm) 
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with hi{t) = ^h{t/V) being a generic distribution with 
zero mean and the variance . The probability Pnm of 
being nonzero for any element is taken to be of the form 
Pnm — if /(I n~m I /&), where f[t) decays exponentially 
or faster at infinity and X^r^o = 1- such a 



definition the mean number of nonzero elements per row 
is just 2M and the quantity S = b/M is natural to term 
" sparsity" . 

In order to calculate the LDOS let us follow the method 
used previously for sparse matrices without band struc- 
ture |2^. Relegating all technical details to a more ex- 
tended publication, let us just mention that the mean lo- 
cal DOS is expressed in the form of a functional integral, 
which is finally calculated by the saddle-point approx- 
imation exploiting the parameters 6^1 and iV ^ 1. 
Unfortunately, it is hard to get rigorous estimate for the 
domain of validity of the method for different values of 
a, M and b. However, for the case of fixed sparsity one 
can show that it is actually exact everywhere in the non- 
perturbative regime of the model. 

Performing this straightforward, but lengthy calcula- 
tion and treating the discrete variable a I""™ in the 
limit 3> 1 as a continuous one, we express the LDOS 
Pl{E — an) in terms of a function G{zi,loi) satisfying 
the equation 



(5) 



where the notations /i — ab/V and zi = z/V are used. 
The LDOS is expressed in terms of the solution G{zi , oji ) 
as follows: 

1 1 f°° 

pUz) - yP{z/Vh p{z,) = -Re / dc^e-^/2-G(.„.) 
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The pair of equations Eqs.(||)-(g) constitute the main an- 
alytical result of the present paper and give the expres- 
sion for the LDOS of sparse WRBM in the closed form. 

Let us first consider the case of the fixed number of 
nonzero matrix elements per row: M — const when 
6 — > GO. At given M the form of LDOS depends on 
the value of parameter /i — ab/V . When this parameter 
is of the order of unity the LDOS form can not be univer- 
sal, i.e. it should be dependent on the particular form of 
the distribution ft.(r) and on M . However, the important 
universality emerges in the limit /i 3> 1. It appears that 
in the large domain of wi , zi the solution of the equation 
Eq.(||) can be approximately written as G{zi,uji) = goWi, 
where go is independent of zi. Indeed, substituting this 
expression in the right-hand side of Eq.(^ we get: 



(7) 



that proves that our approximation is self-consistent as 
long as zi,LUi <C /i. Such a form of G{zi, lui) immediately 
results in the BW form of the LDOS, Eq.(3), in the region 
\z\/V <C fj,, where the width F of the LDOS is given by: 



fi b a 



(8) 
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This allows to give the estimate for the maximal local- 
ization length Imax =T/a (X MV^/a^b. As was pointed 
out above, our nonperturbative treatment (in particular, 
the saddle point evaluation of the functional integral) can 
be valid only if Imax 3> 1. Together with the condition 
fi ^ 1 this gives the following restriction for the region 
of existence of the BW regime: 



< 



'M 



v/2vr/(0)/6 



(9) 



Another case deserving separate investigation is when 
the sparsity S — b/M is fixed at 6 ^ cxd rather then 
M itself. In particular, 5 = 1 case corresponds to the 
standard WRBM where we expect our formulas to be 
in agreement with those derived by Wigner [Q. Since 
& ^ 1 is equivalent for this case to M ^ 1 we can 
search for the solution of equation for G{zi,uji) in terms 
of an expansion with respect to 1/M. One can satisfy 
oneself that at the leading order in 1/M one can put 
G{zi,LUi) — l/2gi^{zi/\/bS)uji\/bS where the function 
gK,{x) is the solution of the following equation: 



9Kix)^--[ f{\u-x\/K) 



iu - g^iu) 



; Reg^{u) < 0. 
(10) 



where the control parameter k = aV ^VbS is intro- 
duced. The LDOS is given correspondingly by 



, , 1 1 

Pl(z) = -Re- „ . 

TT iz — Rgii(z/R) 



i?2 ^ bV^/S (11) 



= ^{iz - 
ll) are identi- 



After introducing the function r{z) 
Rg^{z/R))^^ the two equations Eqs.(|l~ 
cal to those derived by Wigner H] provided a = 1, 5 = 1 
and f{t) is the step function: f{t) = d{\t\ — 1). It is evi- 
dent that in the finite sparsity case the form of LDOS is 
completely determined by the only parameter n oc ab^^'^. 
This fact is in agreement with the mentioned numerical 
observation [|l6|. 

When K <C 1 one can immediately find from Eqs.(|l0|- 
|ll|) that LDOS is given by the standard semicircular law 
Eq.(|l]) with Ro = {SbV^/Sy/'^ as the radius of the semi- 
circle. When K increases the BW form of the LDOS (|) 
emerges , the width F and the domain of validity of BW 
formula being given by the same expressions Eqs.(^^. 
However, now the width F is actually independent of 
the parameter 6 because of b/M = const. One should 
stress that in the BW regime the value of k is bound be- 
cause of the conditio n Eq.(^ which can be rewritten as 
•\/27r/(0)6 ^ K ^ ■\/27r/(0). This domain is quite nar- 
row for not large enough values of b taken in numerical 
experiments |9|JlT[|. 

To show the transition from the semicircle regime to 
that of the BW we have performed the numerical study of 
the equations Eqs.( p^[jTl| ) together with the direct com- 
putation of the LDOS from eigenvectors of large matrices 



of size A'' = 1000 and band size & = 10 without sparsity: 
5 = 1. The parameter V was kept constant in a way en- 
suring the semicircle radius Rq to be unity whereas the 
parameter a used to change k. The function f(\n — m\/b) 
was taken to be a step- function: / = for |n — m| > b. 
For such parameters b^ <^ N and the finite size correc- 
tions are small. On the other hand, still there is a region 
for K where the BW shape for the LDOS is expected. The 
results are presented in Fig.l where the smooth curves 
are solutions of Eqs.( p^-|ll| ) and histograms are obtained 
by the average of Eq.(2p^ovcr number of matrices and 
over those eigenvectors that are not sensitive to the fi- 
nite value of A^. From this figure the whole transition 
from the semicircle to the BW form is seen in detail. 
Unexpected peculiarity of these data is that in the crit- 
ical point (for k w 2) there is a local minimum in the 
center of the LDOS. 

Specific question is about the form of tails of the 
LDOS in the BW regime [|]jl|Jll. It is rather hard task 
to extract the analytical expression for these tails from 
Eqs.( p^[]ll] ). For this reason additional numerical check 
has been done for k = 8 and Rq = 1, see Fig. 2. From this 
figure a sharp transition from the BW to the exponential 
form is clearly seen a.t E — an fa ab — kRo\/8 w 2.8 . 
The detailed analysis shows that it is impossible to dis- 
tinguish between the pure exponential dependence and 
that found in Rcf. where the Wigner's expression for 
the tails [0 was corrected (see also [|l2|). 

The results obtained here allow to understand the na- 
ture of the scaling parameters governing the statistical 
properties of sparse WRBM in different regimes. Indeed, 
statistical properties of eigenvectors and eigenvalues are 
expected to be determined by the ratio (3* = I /Imax 
where I is actual localization length in the unperturbed 
basis (see details in ||l^). Numerical experiments ^-p^] 
for the standard WRBM indeed revealed that the be- 
haviour of the system is completely determined by the 
value of the scaling parameter A = loo/hnax where loo is 
the localization length I for the limit a = 0. This param- 
eter can also be called the "ergodicity parameter" |Tc| ] 
since at A >> 1 any eigenstate is spread uniformly over 
the scale Imax and thus occupies the maximal possible 
number of available states in phase space (in this case 

r(A)«i). 

In it was found that sparsity does not actually af- 
fect the basic fact of proportionality /qo 5^ at a = 0. 
However, the constant Cm — loo/b"^ is expected to be 
very small at small enough values of M. In contrast, 
in the fixed sparsity case M oc & 3> 1 one typically has 
Cm ~ 1. In the latter situation the scaling parameter is 
easily shown to be A oc (a5)^ for the BW regime. Then 
the condition k ^ 1 is equivalent to A ^ 6, immediately 
showing that the BW regime corresponds to the com- 
plete ergodicity: loo ^ Imax-^ and as a consequence to 
the Wigner-Dyson level statistics (see discussion in [p2[). 
However, if the number M is fixed rather than the ra- 
tio b/M, one finds A oc CM(a5^/^)^ for the BW domain 
restricted by A Cm&- For very sparse matrices the 
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latter condition is compatible with " non-ergodicity" cri- 
teria A ^ 1 in view of Cm ^ 1- Indeed, in the paper 
where the extreme case M = 1 was studied, the level re- 
pulsion was found to be quite weak and dependent on the 
combination ab^^^ in agreement with the argumentation 
presented above. 
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be generic for the RBM with diagonal elements having 
zero mean value but strongly fluctuating in comparison 
with the off-diagonal ones. The detailed discussion can be 
found in: Ph.Jacquod, D.Shepelyansky Phys. Rev. Lett., 
v. 75 (1995) 3501; Y.V. Fyodorov, A.D.Mirlin , Phys.Rev 
B v.52 (1995) R11580; K.Frahm, A.Miiller-Groeling , Eu- 
roph.Lett. v.32 (1995) 385 
[21] A.D.Mirlin and Y.V.Fyodorov, J.Phys. A 24 (1991) 2273 
[22] Y.V.Fyodorov, A.D.Mirlin, H.-J.Sommers, J.Phys. 
(France) 2 (1992) 1571. 

FIGURE CAPTIONS 

Fig.l Dependence of the LDOS on the control param- 
eter K. 

Fig.2 The tail of the LDOS for k = 8 and Rq = 1. 
The smooth solid curve is the solution of Eqs.([lC|-pTl); 
the dotted curve (a) is the Breit-Wigncr formula (^; the 
dotted straight line (b) is the exponential dependence. 
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